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ABSTRACT 

The importance of magnetic fields in three-dimensional magnetoconvection models of the Sun’s 
photosphere is investigated in terms of their influence on the continuum intensity at different viewing 
inclination angles, and on the intensity profile of two [O i] spectral lines. We use the RH numerical 
radiative transfer code to perform a posteriori spectral synthesis on the same time-series of magneto¬ 
convection models used in our publications on the effect of magnetic fields on abundance determination. 
We obtain a good match of the synthetic disc-centre continuum intensity to the absolute continuum 
values from the FTS observational spectrum; the match of the centre-to-limb variation (CLV) syn¬ 
thetic data to observations is also good, thanks, in part, to the 3D radiation transfer capabilities of the 
RH code. The different levels of magnetic flux in the numerical time-series do not modify the quality 
of the match. Concerning the targetted [O i] spectral lines, we find, instead, that magnetic fields 
lead to non-negligible changes in the synthetic spectrum, with larger average magnetic flux causing 
the line to become noticeably weaker. The photospheric oxygen abundance that one would derive 
if instead using non-magnetic numerical models would thus be lower by a few to several centidexes. 
The inclusion of magnetic fields is confirmed to be important for improving the current modelling 
of the Sun, here in particular in terms of spectral line formation and of deriving consistent chemical 
abundances. These results may shed further light on the still controversial issue regarding the precise 
value of the solar oxygen abundance. 

Subject headings: magnetohydrodynamics (MHD) — radiative transfer — Sun: abundances — Sun: 
granulation -- Sun: magnetic fields — Sun: photosphere 
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1. INTRODUCTION 

Sun is magnetised over 
even in the so called 


its whole sur- 
' q uiet” r egions 


Trujillo Bueno et al. 200 1 : ISolankil 


Sanchez Almeida & Martinez G onzalez! _ 1 _ - __ 

Bellot Rubio fc Orozco Suared 120121 iNordlund et al l 


200^), with at least as much magnetic flu x present in 
the quiet S un itself as in act ive regions dGosic et all 
l2014t [Thornton fc Parnellll2011l) . The presence of ubiq¬ 
uitous magnetic fields has an influence on the average 
temperature stratification in the solar photosphere, 
thus indire ctly affecting the many spe ctral lines formed 
there (see iFaLbian et all 120101 I2012L called Papers I 
and II, respectively, in the following). This, in turn, 
can have an impact on the accuracy of the chemical 
composition derived from different spectral features, 
when non-magnetic atmospheric models, even current 


geneity (e.g., Stein fc Nordlund 1998c As 

plund et al. 

2000 

: iFrevtag et al. 20021 Asplund 2005| 

Scott et al. 


2014), are employed for solar abundance derivation 


purposes. 

Three-dimensional (3D) convection simulations have 
already been extensively applied to aid in the under¬ 
standing of different solar phenomena, including the im- 
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pact of the convective velocity field on the formation of 
spectral lines in terms o f their widths , shift s, shapes 
and a symmetries (e. g.. IDravins et all 119811: IDravinsI 
Il982t INordlund et all 20091) . Generally, though, the 
effects of magnetic fields on the determination of so¬ 
lar chemical abundances has not been systematically 
explored. In Papers [| and [U, we carried out a first 
quantitative assessment of the impact of magnetic fields 
in 3D convection models on the determination of the 
solar chemical composition. To that end, we calcu¬ 
lated time-series of magnetoconvection simulations us¬ 
ing the Copenhagen STAGGER radiation magnetohy¬ 
drodynamics code. The quality of those simulations 
was confirmed in follow-up studies, where the snapshots 
were te sted for thermodynamic al and polarisation prop¬ 
erties dBeck et al.l 120131 120151) . as well as being em¬ 
ployed for associating G-band bright po ints with mag¬ 
netic field strengths of a few kG dCriscuoli fe UitenbroeU 
l20l . and for comparing numerically-predicted and 
observationally-derived velocity fields in sol ar convec¬ 
tion via Fourier power spectra (lYelles Chaouche et all 
120141) . Our 3D magnetohydrodynamical (MHD) mod¬ 
els were also used for the interpretation of solar irra- 
diance measurements, whereby the predictions of the 
simulations may provide a possible explanation for the 
counter-phase variation of the spectral irradiance at vis- 
ibl e and IR wavelengths with the solar activity cycle 
(jCriscuoli fc Uitenbroekll2014l) . 










































































2 


The motivation for the present paper is to further 
test the quality of our magnetoconvection models and 
to proceed in the investigation of the effect of magnetic 
fields on the average spectra and derived chemical 
composition of the solar photosphere. Using the MHD 
simulations of Papers U and [I]], 

(1) we attempt an absolute match of the synthetic con¬ 
tinuum radiation to disc-centre FTS observations; 

(2) we carry out a study of the continuum intensity vari¬ 
ation when moving from disc centre to limb; and 

(3) we investigate whether those models also correctly re¬ 
produce selected [O i] absorption lines in the solar spec¬ 
trum. 

Concerning item (1), iTruiillo Bueno &; Shchukina 
(2009) used a single 3D HD snapshot by lAsplund et ah 
1(2000( 1 to produce a synthetic disc-centre continuum in¬ 
tensity spectrum. Values at wavelengths in the visible 
part of the spectrum proved to be roughly 10% below 
observations. In Papers Q] and ITU, using 1.5D radiation 
transfer calculations, we obtained a match of the abso¬ 
lute values of the synthetic continuum intensity at disc 
centre to within ~ 7% of the solar observational data. 
In the present work, we see how this match can be sub¬ 
stantially impr oved by using a d i fferen t spectral synthe¬ 
sis code (H I I. iuitenbroekll200ll [2002D . with, therefore, 
opacities different to those used in our previous papers. 

Conce rning item (2), i . e. the continuum CLV be¬ 
haviour, Koesterke et al.l (1200811 found tha t shortcomings 
in CLV results from the 3D HD model of lAsplund et al.l 
(1200011 were roughly comparable to the shortcomings of 
classical ID models. The results indicated the pres¬ 
ence of too steep a temperature gradient in the 3D HD 
model. They tentatively concluded that the improved 
description of the convective energy transport in 3D 
may be offset by deficiencies introduced by a too sim¬ 
plified radiation transfer scheme. Including more opac- 
ity bins and additio nal sources of continuum opacities led 
IPereira et al.l (|2013l l (called Paper [Hi] in the following) to 
achieve very good agreement of their CLV results with 
observations. It is thus of interest to check how well our 
3D MHD models match the observations concerning the 
continuum CLV behaviour, in spite of the reduced num¬ 
ber of opacity bins, by using the 3D radiation transfer 
capabilities of the RH code. 

In the final part of the current work, we will look at the 
effect of the magnetic fields in convection models on the 
formation of selected [O i] spectral lines. IPereira et al.l 
have recently argued in Paper [Hl| that claims for revised 
abundances based on 3D MHD mod elling are premature. 
On the other hand, iSocas-Navarrol (12014T) very recently 
found an extreme sensitivity of the derived solar oxy¬ 
gen abundance to possible systematics (retrieving val¬ 
ues in the range loge(O)© = 8.70 — 8.90 dex). Typical 
uncertainties, mainly intrinsic to employing atmospheric 
models, would thus be substantially higher than what 
is usually claimed in the corresponding literature. It is 
clear that the solar oxygen chemical abundance deriva¬ 
tion is particularly difficult, because the corresponding 
absorption lines present in the solar spectru m are few 
and problematic (see review bv lAsplundll2005[ ). 

In this work, we will show that, barred the remaining 
uncertainties (on, e. g., continuum opacity, treatment of 
radiative transfer and non-local thermodynamic equilib¬ 


rium effects) that complicate the treatment of spectrum 
formation, as well as t hose a ffectin g the solar spectrum 
observations (see, e. g.. iBeck et alJI20lD h the MHD ap¬ 
proach provides synthetic spectra that match very well 
both the solar continuum intensity at different inclination 
angles, as well as the [O i] 557 nrn and 630 nm absorption 
lines. 

In the following section, we explain the simulation 
setup, we inform on the input data chosen for the spectral 
synthesis computations, and briefly present the method 
adopted in this study. In sections mm and[5j we present 
the analysis of the theoretical results of interest on disc- 
centre absolute continuum intensity, on centre-to-limb 
variation of the continuum intensity, and on continuum 
intensity normalised spectral line profiles, respectively, 
and compare them to available observational data. The 
paper ends with a discussion and conclusions (Sec. [6] and 
0 . 

2. SIMULATION SETUP AND RESULTING MAGNETIC 
FIELD STRATIFICATION 

2.1. Setup of the numerical (M)HD simulations 

The 3D model atmospheres adopted as input for the 
spectral synthesis are based on snapshots selected from 
the statistically stationary regime of the convection in 
our previous (M)HD simulations (PapersQjandQf]), which 
used input solar abundances of carbon, oxygen and iron 
of, respectively, loge(C) Q = 8.55, loge(O)© = 8.77, and 
loge(Fe)© = 7.50. The snapshots we employ in this 
study belong to the different MHD series in those papers, 
which had a vertical magnetic field of uniform strength 
implanted at time t = 0 of 56 G, 112 G and 224 G, re¬ 
spectively (labelled “50 G”, “100 G” and “200 G” in the 
following). Additionally, we will also use here the “0 G”, 
i.e., purely HD , series of those papers. 

The atomic and molecular transition data we em¬ 
ploy in the intensity spectrum calculations are from the 
Vienna Atomic Line Database ( VALD) and from the 
National Institute of Standards and Technology Atomic 
Spectra Database ( NIST ASP). VALD and its deve l¬ 
opment have been descri bed b v iPiskunov et ahl (11995D . 
iRvabchikova et all (|1999t ). and iKunka et alJ (1 1 9991) . In 
particular, we made use of the recent (February 2014) 
third rele ase of VALEB, and of NIST ASE@ in its online 
version 5 (jKramida et al.ll2012D . 

We empl oyed version 2 of the RH code (jlJitenbroekl 
I2001L120021) to perform the 3D radiation transfer calcu¬ 
lations (we refer to them as “R3D” in the rest of this 
article) on the basis of the input 3D (M)HD snapshots. 
In the solar atmosphere, both Rayleigh scattering (off 
neutral hydrogen and helium) and Thomson scattering 
(off free electrons) may cause a departure from LTE con¬ 
ditions even in the absence of NLTE spectral lines. These 
scattering effects are negligible for the continuum at disc 
centre, but may become progressively more important 
towards the limb, where higher layers are being sam¬ 
pled. We thus include Rayleigh-Thomson (RT) scatter¬ 
ing in the calculations of this paper, with a view to the 
centre-to-limb variation calculations of Sect. 0 Verti¬ 
cally, we cut the input model to 64 depth points covering 
a physical depth of -0.5 to 0.5 Mm approximately. Hori- 

2 http://vald.astro.univie.ac.at/ vald3/php/vald.php 

3 http://physics.nist.gov/asd 
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zontally, the full sampling of the original snapshots was 
maintained for the centre-to-limb variation (CLV) runs, 
with 252 grid points over 6 Mm in each direction, while 
the sampling was reduced to 63 points in each horizontal 
direction over the same physical extent for the spectral 
line synthesis calculations. 

We first computed the absolute intensity spectrum at 
disc-centre. Then, we performed the CLV runs using a 
Gauss-Legendre set with 9 inclination angles and 4 az¬ 
imuth angles per octant (thus providing 36 angles per 
octant, or 144 rays in total). In sections [3] and 0] we will 
be comparing the RH results of this paper with those 
of simpler (e. g., 1.5D, or without RT scattering) radi¬ 
ation transfer calculations carried out with other codes: 
the corresponding details will be explained in those sec¬ 
tions. For the spectral line calculations, which focused 
on the 557 nrn and the 630 nrn oxygen features, the ra¬ 
diative transfer problem was solved using a fine (10 mA) 
and uniform sampling in A space. The resulting output 
spectrum was set to cover 1.2 A around each of the two 
oxygen spectral lines of interest (which allowed us to ac¬ 
count for any significant blends affecting them). 

The spectral synthesis results for the MHD case were 
obtained neglecting the spectral line broadening due 
to the Zeeman effect, on the basis of the expectation 
that, compared to the indirect effect due to temperature 
changes, this direct effect of magnetic fields on spectral 
lines is relatively small. It in fact increases with A 2 and it 
thus generally becomes appreciable in the infrared (IR), 
rather than in the visible wavelength range where the 
spectral lines we selected for this study fall. 

2.2. Average magnetic field distribution in the 
convection cells 

Some of the results presented in this paper are associ¬ 
ated with the introduction of magnetic field in convec¬ 
tion models, which can induce significant stratification 
changes, with consequent effects on spectral lines. Thus, 
for a better appreciation of our findings, it is adequate to 
consider briefly the average values of the magnetic field 
in the snapshots that serve as a basis for the radiation 
transfer calculations. The upper panel of Fig. [I] shows, 
using solid curves, the height distribution of the average 
magnetic field intensity, with the average being carried 
out separately for each horizontal plane and along time; 
the figure, therefore, shows the discrete version of: 

dt J dx dy B(x,y, z,t) , 

horizontal 
plane 

(!) 

with A the area of the horizontal plane in the numerical 
box and to and t/ the limits of the time interval used 
for the averaging. Here, and for all subsequent figures in 
this subsection, the time averaging is done using the last 
100 snapshots (corresponding to 50 solar minutes) of the 
statistically stationary regime in the simulations, with 
snapshots separated by 30 solar seconds. The different 
series are identified with blue, green and red colour for 
the “50 G”, “100 G” and “200 G” series, respectively. All 
three ( B) distributions decrease smoothly with height in 
the topmost Megameter of the convection zone and reach 
values of ~ 80 G, ~ 150 G, and 260 G, respectively, at 


(B)(z) = 


1 1 
tf — to A , 


to 


z = 0. In the highest layers (i.e., roughly in the top¬ 
most 150 km), the distributions show an upturn, with 
the value of the mean magnetic held strength increas¬ 
ing to a maximum of ~ 100 G, ~ 170 G, and ~ 300 G, 
respectively for the “50 G”, “100 G” and “200 G” series. 

It is worth noting that the “100 G” (H)-distribution 
of Fig. [1] is_ n ot fa r from the di stribution shown 
by iShchukina fc Truiillo-Buenol (|2Qlll . i2015f) to match 
the scattering polarization signal (induced by the 
Hanle e ffect) from a heig ht of around 250 — 300 km 
(jGurtovenko fc Kostvkl fl989h . as derived from observa¬ 
tions of the Sri 460.7 nm line, simultaneously allowing 
for the explanation of the polarization signal (induced by 
the Zeeman effect) from a height of around 60 km, as de¬ 
rived from observations of the Fel 630.25 nm li ne. T he 
small-scale local dynamo models of iVogler fc Schiisslerl 
lt200l or iRempel (2014j ). instead, lead to too weak av¬ 
erage magnetic held strength values to explain the Sr I 
460.7 nm Hanle-effect data. 

The upper panel of Fig. [T] also shows (thick dashed 
curves) the height dependence of the root mean square 
(rms) magnetic held intensity, B rms = yj(B 2 ), which 
also serves to illustrate the distribution of magnetic en¬ 
ergy in the box in the stationary regime. The B rms 
curves peak (with values of ~ 300 G, ~ 460 G, and 
~ 670 G, respectively) in progressively higher layers (be¬ 
tween ~ 1700 and ~ 1200 km) as magnetic hux increases 
between the “50 G” and the “200 G” case, then decrease 
monotonically. 

The dotted horizontal lines in the hgure mark the con¬ 
stant value of the vertical component of the held that was 
implanted at t = 0 in each series. In all cases, in the sta¬ 
tionary regime the B rms held is larger than the average 
held intensity. This is a result of the intermittency cre¬ 
ated in the held by the action of the convective hows. To 
show this effect from a different perspective, we plot in 
the lower panel of Fig.[l]the ratio yj ( B 2 )/ ( B). This ratio 
can serve as an index for the intermittency level in each 
of the time-series. The ratio is seen to increase monoton¬ 
ically with height in the subsurface layers, and to reach a 
maximum very near the surface, z = 0. This is expected 
since the Alfven Mach number of the convection hows 
is largest near the surface. We also see that the inter¬ 
mittency level decreases the larger the initially implanted 
held, with peak ratios of ~ 3.0, 2.4 and 2.1 for the “50 G”, 
“100 G” and “200 G” series, respectively. Above the sur¬ 
face, the intermittency ratio decreases sharply, given the 
less passive role of the magnetic held generally for higher 
atmospheric levels and the decreasing intensity of the 
convection there. 

3. DISC-CENTRE ABSOLUTE CONTINUUM INTENSITY 

In this section we present the hrst group of results in 
the paper, namely the comparison with observations of 
the absolute disc-centre continuum intensity values ob¬ 
tained through spectral synthesis calculations from our 
(M)HD snapshots (Fig. [21 upper panel). For the observa¬ 
tions we took the Fourier Transform Spectrograph (FTS) 
data, drawn in the hgure as a collection of points in light 
gray. For the comparison, a pseudocontinuum is needed, 
to use as an envelope to the FTS data. We obtained 
it (black line) by linking the local maxima of the orig¬ 
inal FTS data, bridging gaps across wavelength regions 
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Figure 1 . Upper panel: height-stratification of the average magnetic field intensity (solid curve) and of the root-mean-square magnetic 
field (thick dashed curve), expressed in Gauss units. Blue colour is used for the “50 G” series, green for the “100 G” series, and red for the 
“200 G” series. The dotted horizontal lines mark the constant value of the horizontally-averaged vertical component of the signed magnetic 
field intensity for the three series. Lower panel: height-stratification of the ratio of the rms magnetic field intensity, B rrns , and of the 
average magnetic field intensity < B >. Same colour coding as for the upper panel. 


affected by strong absorption features. Our R3D data 
are shown as blue diamonds and they represent the abso¬ 
lute continuum intensity values obtained using the three- 
dimensional radiation transfer calculation explained in 
Subsection O To take into account the fluctuations in 
time, we calculated the value of the continuum at each of 
the wavelengths as the time average of 21 snapshots cov¬ 
ering 50 solar minutes, with the first snapshot taken at 
100 solar minutes inside the statistically stationary part 
of the simulation. The theoretical values (the blue dia¬ 
monds) are very near the FTS continuum, as also appar¬ 
ent in the lower panel of the figure: the average mutual 
deviation is ~ 3.8%, with the minimum and maximum 
relative difference being ~ 1.5% (at 600 nm) and ~ 6% 
(at 900 nm), respectively. 

In the upper panel of Fig. [2j we also show, as red tri¬ 
angles, our previous results of Paper [H, obtained with 
a different radia tion transfer ca l culati on, namely, using 
the Lilia code of iSocas-Navarrd (|2001il and without any 
RT scattering. The better match to observations of our 
current results compared with those of Paper [jl] provides 
reassurance both on the quality of our models and on the 
opacities used by the RH code. We further tested the ef¬ 
fect of not including Rayleigh and Thomson scattering 
in our current RH calculations. We found the change in 
absolute continuum intensity to be unimportant. 

The small remaining percentage difference (shown in 
Fig. m lower panel) of the R3D results with respect 


to FTS can be understood as partly due to the fact 
that the (M)HD snapshots have a slightly lower (within 
< 20 K) outgoing radiation flux compared to the mea¬ 
sured solar one. However, that can explain only < 1% 
of the discrepancy. A more significant contribution to 
the ~ 3.8% average difference with FTS may arise from 
the comparatively low number of opacity bins used in 
our Stagger-code calculation. A thorough treatment of 
line blocking/blanketing effects by using very accurate 
atomic data for many spectral lines may also provide a 
further improvement of the match with observations, an 
approach which can, however, become prohibitively time- 
consuming when coupled with 3D radiation transfer, due 
to the plethora of lines present in the solar spectrum that 
need to be taken into account. In any case, the FTS atlas, 
like any observation, is, itself, affected by uncertainties 
of up to a few percent (for example af fecting r egions 
of str ong telluric absorption, see IBurlov-Vasiliiev et al.l 
(119981V ). which contribute to part of the difference seen 
here. 

The comparison between the theoretical results for HD 
and MHD runs is also of interest. It turns out that there 
is only a small difference between them; in fact, the HD 
curve would be nearly indistinguishable if plotted along¬ 
side the blue diamonds of Fig. [2] (which correspond to the 
R3D MHD results for the “200 G” case): the HD contin¬ 
uum is only A (I) < 0.02 • 10 4 W m~ 2 sr nm below the 
MHD one. The reason why their continuum level turns 
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A [nm] 

Figure 2. Absolute disc-centre continuum intensity at different wavelengths (upper panel). The solar atlas data from FTS observations 
are shown in light grey colour. The theoretical R3D values are plotted as blue diamond symbols. Our FTS pseudocontinuum is shown as a 
solid black curve. Our R1.5D values of Paper [nj are shown as red triangles, fn the lower panel, we show the percentage difference between 
the R3D results and the FTS pseudocontinuum. The dashed horizontal line marks the ~ —3.8% average (R3D MffD - FTS) deviation 
level. 
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out to be nearly identical is that (Paper [J, Fig- 3; Paper 
in, Fig. 1) the MHD-HD temperature difference tends to 
be smallest (tending to negligible) precisely in the layers 
where the continuum forms. 


4. CENTRE-TO-LIMB VARIATION OF THE CONTINUUM 
INTENSITY 

As an observer’s line-of-sight moves towards the solar 
limb, the continuum intensity decreases because it orig¬ 
inates in higher layers which are not as bright (warm) 
as the deeper ones to which the observer’s line-o f-sig ht 
penetrates when looking at disc centre (e.g.. !Gravlll992l 
chapter 9). We can thus use the behaviour of the theoret¬ 
ical CLV to verify whether the temperature stratification 
of the models results in a proper matching of the solar 
continuum intensity CLV. 

In Fig. [31 we show our results on continuum intensity 
centre-to-limb variation (CLV) for the synthetic spectra 
normalised to the corresponding value at disc centre. We 
studied the continuum intensity behaviour both against 
p, z = cos$hciio (upper panel), with (?heiio the heliocentric 
angle, as well as against wavelength (lower panel) in the 
range 600 nm to 2400 nm. In the upper panel of Fig. [31 
the theoretical R3D continuum intensity results from our 
computations for the MHD case are represented as dia¬ 
mond symbols linked by dashed s egments, while obser¬ 
vation^ are shown as ’+’ signs for iNeckel fe Labs! (1199411 
and ’x’ symbols, for iPierce et all ( 1977f) . Our R3D re¬ 
sults are seen to match very well the solar limb darkening 
observational data. Taking into account radiation trans¬ 
fer effects properly in a 3D inhomogeneous atmosphere 
thus produces realistic values of the average radiation in 
the continuum. In the figure we are also showing (as 
triangles linked by a dotted line) the results obtained us¬ 
ing a simplified radiation tran s fer c alculation with the 
Nicole code of iSocas-Navarro et all ( 2014 ) (which, like 
the Lilia code before, uses a different opacity package 
than RH). The Nicole version used had a 1.5D approach 
to the horizontal radiation transfer (called “R1.5D” in the 
following). In that approach, we considered as many p 
angles as in the full 3D calculation just shown, but the 
horizontal radiation transfer coming out of a given pixel 
is calculated via an infinite repetition of the correspond¬ 
ing column in the two horizontal directions. The R1.5D 
results are found to deviate substantially from the ob¬ 
servational data: differently from the previous Section’s 
result on disc-centre continuum intensity, the continuum 
CLV values are not only sensitive to adopted opacities, 
but are also sensitive to horizontal transfer effects, espe¬ 
cially as one moves toward the limb. It is then natural 
that the R1.5D results deviate importantly from the ob¬ 
servations when approaching the limb, since they do not 
take proper account of the hills and valleys of the r iso¬ 
surfaces in the 3D models. 

The lower panel of Fig. [31 compares the observational 
data and theoretical results just discussed, now against 
wavelength. This more clearly highlights the differences 


4 Note that the INeckel fc Labd (119941 ) data agree well with the 
less recent ones by IPierceTT*5!aughteii (1 197711 (not shown here) 
in the wavelength interval of 303 — 730 nm common to the two 
datasets. In the adjacent c ommon wavelength interval of 740 — 
1100 nm, th e Neckel & Labs (1994) data agree well with those of 
I Pierce et all <119771 1. For this reason, we show here the values for 
the latter dataset only for the non-overlapping wavelength region 
above 1100 nm. 


between the R3D and R1.5D approaches. The symbols 
are the same as in the upper panel of this figure (dia¬ 
monds for R3D, triangles for R1.5D, ’+’ signs and ’x’ 
for the observations), except without colour coding in 
wavelength. Again here, the R3D-R1.5D deviations (i. e., 
those due to horizontal transfer effects) generally increase 
when moving toward the limb, in particular, at longer 
wavelengths. This is in line with the results in Paper 
153 (their Fig. 3), where a similar behaviour is visible in 
the case of the classical MARCS and PHOENIX ID so¬ 
lar photosphere models (while the ID Holweger & Muller 
model, being semiempirical, not surprisingly matches the 
observational data as well as their 3D models). It should 
be noted that the comparison of theoretical and obser¬ 
vational data very close to the limb (i.e., at p z < 0.2 or, 
equivalently, 0} le i; o > 80 degrees) is less helpful, since the 
latter data tend to become less reliable there. 

In both panels of Fig. [31 we show only the “200 G” 
MHD case, for both the R.1.5D and R3D radiative trans¬ 
fer calculations. The R3D case for the other series, in 
fact, generally displays a very similar match to observa¬ 
tions as the one plotted. The situation is thus similar to 
what was already seen in the previous Section regarding 
the absolute continuum intensity, with relatively small 
differences between the results for the HD case and the 
MHD case(s). 

The fine match to observations of our current R3D re¬ 
sults again provides reassurance on the accuracy of our 
3D (M)HD models and on the RH opacities and spectral 
synthesis. Differently from Sec. [31 though, it is harder 
to disentagle for these CLV runs how much of the differ¬ 
ence between R3D and R1.5D results is due to, respec¬ 
tively, different opacities, Rayleigh and Thomson scatter¬ 
ing (which should play more of a role at inclined angles), 
and pure geometrical / horizontal transfer (3D vs. 1.5D) 
effects in the treatment of radiation. 

5. [O I] SPECTRAL LINE PROFILES 

In this section, we analyse the behaviour of the pro¬ 
file of selected spectral absorption lines when using 3D 
snapshots from our HD and MHD series and when chang¬ 
ing the C and O solar abundance adopted as input to 
the spectral synthesis code. The main aim is to under¬ 
stand how the change in temperature stratification due 
to the presence of magnetic fields affects the spectral lines 
and to additionally verify whether our 3D MHD photo- 
spheric models coupled with 3D radiative transfer for the 
a posteriori spectral synthesis can fit observations bet¬ 
ter, or at least as well as when using input photospheric 
models from our HD series. For this part of our study, 
we selected two representative spectral lines of neutral 
atomic oxygen, namely [O i] 557 and [O i] 630 nm, both 
produced by radiatively forbidden transitions. These 
two spectral lines are deep-form ing and very weak in 
the_solar spectrum. According to iGurtovenko fe Kostvkl 
(11989D . the main contribution to the formation of their 
core is from a range of layers around a height of, re¬ 
spectively, ~ 130 km and ~ 180 km, above 7500 nm = 1. 
Being formed in LTE, and given the scarcity of avail¬ 
able oxygen lines in the solar spectrum, these spectral 
lines represent important tools in terms of t he solar 
oxygen abundance dete rmination (Melendez fc Asplundl 
120081 : [Scott et all 120091 ). Their use for the purpose of 
estimating the solar oxygen abundance is, however, sig- 
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Figure 3. Upper panel: predicted continuum intensity at different wavelengths (indicated by the different colours, as labelled) from our 
R3D (this paper, diamond symbols linked by dashed segments) and R1.5D (Paper II, triangular symbols) comp utations f or diffe ren t values 
of fi z = cos 6 (w here 6 is the helioc entric angle) is compared to solar limb-darkening observational data by INeckel Sz Labs! d 19941) (’+’ 
symbols), and bv IPierce et al.l d 19771) (’x 5 signs). The dotted horizontal line indicates the normalisation level with disc-centre continuum 
intensity. Lower panel: the same, now versus wavelength and with the different curves being for different fi z = cos 6 values, as labelled. 
Symbols are the same as for the upper panel, as marked in the legend. 
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Figure 4. Solar observational disc-centre intensity data from FTS (black open circles) and R3D synthetic intensity profiles (coloured 
curves) in HD and MHD for the blended feature at 630 nm, which contains the [O i] line of interest. As per the legend in the figure, the HD 
and MHD synthetic profiles shown are the result of adopting different choices for the input C and O solar abundances. The wavelength of 
the [Oi] line and of the blending Ni i line are marked with a dashed vertical segment and with a solid vertical segment, respectively. The 
dotted horizontal line indicates the normalisation level with the adjacent continuum intensity. 


nificantly affected by blends, mainly C 2 molecular lines 
for the [Oi] 557nm absorption line, and a Nil atomic 
line for [Oi] 630nm. Additionally, the blending Nil line 
requires proper treatment of its isotopic splitting, ac¬ 
counting for which significantly improves the fit to ob¬ 
servational data. In the following analysis, we include all 
significant blends, among them the C 2 and Nil lines just 
mentioned. Note that, while free of significant known 
blends, the strong, widely used neutral oxygen triplet 
lines at 777 nm are prone to significant non-L TE effects 
(see iFabbian et al.1120091 iPereira et al.1 [2009) and thus 
would require a more complicated and time-consuming 
line synthesis treatment than the one we adopt here. 

The [0 1 ] 557 and [0 1 ] 630 nm lines fall in the wave¬ 
length region where the absolute continuum intensity 
from our theoretical R3D calculations best matches the 
FTS solar atlas (see Fig. 0. We here compare the 
continuum-normalised theoretical spectra with Fourier 
Transform Spectrograph (FTS) observations of the solar 
spectrum. We normalised the spectral intensity across 
each of the wavelength regions to the relevant intensity 
of the adjacent pseudocontinuum. Note that for spectral 
lines around 600 nm, the variation of the continuum level 
with wavelength across the line is negligible (< 10~ 4 ) and 
does not need to be taken into account for our purposes. 

The line formation calculations are computationally 
demanding, since many wavelength points must be used 
(as explained in Sec. O) for an accurate study of the 
line profiles. Therefore, we selected as input model pho¬ 
tosphere a single snapshot per series, namely one from 


the series with no magnetic field and another one from 
the “200 G” MHD series, respectively labelled “HD” and 
“MHD” in the following. Both snapshots were taken at a 
point in the calculations which corresponds to 100 min¬ 
utes of solar time evolution in the statistically station¬ 
ary regime of the convection. The choice to study the 
relevant effects based on the “200 G” MHD series was 
made to estimate the expected maximum changes that 
the presence of magnetic fields induces - via a modi¬ 
fied atmospheric temperature structure - on the profile 
of the selected oxygen spectral lines. Given the weak¬ 
ness of the latter, this also makes the representation of 
results in the figures of this Section clearer. This study 
thus covers quiet Sun to weak plage conditions. 

In Figs. 2] and [5] we show the solar disc-centre inten¬ 
sity results of the R3D line synthesis for, respectively, 
the spectral region around the [0 1 ] 630 nm and 557 nm 
absorption features. Our initial approach was to at¬ 
tempt to fit the FTS observations with the R3D synthetic 
spectra obtain ed using input abundance values from 
lAsplund et al.1 (120091 ). in particular, log e(C) Q = 8.43 dex 
and loge(0)Q = 8.69 dex. The solar carbon abundance 
adopted is important here, since it affects the right calcu¬ 
lation of the C 2 concentration, which in turn is relevant 
due to the corresponding molecular blends affecting the 
[Oi] 557nm spectral feature. The formation of CO, in¬ 
stead, takes place at cooler temperatures than relevant 
here and is therefore unimportant for our results. 

By comparing the corresponding HD and MHD syn¬ 
thetic intensity profiles (red curves in Figs. [Hand [5] re- 
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[0!] 557.734 nm 



Figure 5. Solar observational disc-centre intensity data from FTS (black open circles) and R3D synthetic intensity profiles (coloured 
curves) in HD and MHD for the blended feature at 557 nm, which contains the [O i] line of interest. As per the legend in the figure, the HD 
and MHD synthetic profiles shown are the result of adopting different choices for the input C and O solar abundances. The wavelengths 
of the [O i] line is marked with a dashed vertical segment, and the wavelengths of important blending spectral lines are marked with solid 
vertical segments (C 2 (l) = C 2 557.733 nm; C 2 (2) = C 2 557.735 nm; C 2 (3) = C 2 557.740 nm; Yi = Yi 557.742 nm). The dotted horizontal 
line indicates the normalisation level with the adjacent continuum intensity. 


spectively for [Oi] 630 nm, and [Ol] 557nm), one clearly 
sees that both of the spectral absorption features we tar¬ 
geted become weaker in MHD, as expected due to the 
higher average temperature in the relevant line-forming 
layers. On the other hand, for both spectral lines, the 
fit against observations is not satisfactory when using 
those abundances. The lines are significantly too weak, 
in both MHD and HD. Not including blends would of 
course produce spectral profiles having the wrong shape 
because they would miss the additional asymmetry. 

The poor fit of the red curves in Figs. 0] and [5] is 
to be expected because our (M)HD models have a differ- 
en t tem perature stratification from the ones employed by 
[Asplund et all ( 2009 ). One may thus change the value 
of the elemental abundances adopted as input in the 
spectral synthesis, to verify the sensitivity of spectral 
features to a different chemical composition. We tested 
adopting different C and O input solar abundances and 
visually obtained a best fit when choosing the values 
loge(C) 0 = 8.50 dex and loge(O) 0 = 8.73 dex. The 
green dashed curves in Fig. [4] and Fig. [5] show the pro¬ 
file of the targetted spectral lines when adopting these 
as input abundance values. The thick dashed curve rep¬ 
resents the MHD case, while the thin curve with crosses 
is for the HD case. As seen in the figures, for both of 
the wavelength regions studied here, our synthetic lines 
for the MHD case match very well the disc-centre inten¬ 
sity profiles observed from the Sun, with the fit being 
particularly good in the case of the 630 nm spectral ab¬ 


sorption feature. It is noteworthy that this good agree¬ 
ment in both wavelength regions is achieved while using 
the same set of C and O solar abundance values, as one 
should require for consistency. Note that the ab unda nce 
of Ni w as maintained fixed at the value of lAsplund et al.1 
(|2009f) . loge(Ni) 0 = 6.22 dex. 

In the case of the 630 nm absorption feature, a cru¬ 
cial contribution to the composite intensity profile comes 
from accounting for the isotopic splitting of a blending 
neutral nickel feature present on its longer-wavelength 
side. We achieved this by scaling the Johansson et al. 
(2003) measured log gf for the nickel spectral line un¬ 
der consideration by the abundance of each of five_Nii 
isot opes (the same treatment adopted by iPereira et al.1 
120091 see their Sect. 4.2.2 and their Table 4). Regarding 
the 557 nm absorption feature, the presence of a small 
residual discrepancy - with respect to the FTS data - on 
the red side (at ~ 557.738 nm) of our R3D line profile 
obtained using our final adopted abundances, and seen 
in Fig.0 may indicate that a weak, yet unknown blend 
might be present at that w aveleng th, just redward of 
the [O i] line. We note that iMelendez fc Asplundl (I2008T I 
obtained a very good fit of the 557 nm wavelength re¬ 
gion, at least in the case of the Delbouille solar atlas, 
when using their own spectrum synthesis code on a three- 
dimensional hydrodynamical model of the solar atmo¬ 
sphere obtained by temporally a nd sp atially av eragin g 
the snapshots from the series of lAsplund et al.1 ( 2004 1. 
They derived values of logec ~ 8.43 based on fitting 
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the blending C 2 lines in the relevant wavelength region. 
With that value, they then found logeo ~ 8.70 from 
the 557 nm absorption feature. Those abundance values 
are significantly (~ 0.03 — 0.07 dex) lower than our final 
adopted input C and O solar abundances. 

6. DISCUSSION 

In this study, we have obtained a good match of quan¬ 
tities resulting from our 3D (M)HD simulations against 
observational data. The studied quantities were: the ab¬ 
solute continuum intensity at disc centre, the CLV of 
the continuum intensity, and the profiles for two selected 
absorption features, each containing an oyxgen spectral 
line ([O 1 ] 557 nm and 630 nm, respectively) blended with 
lines of other chemical elements. 

The continuum intensity (disc-centre and/or CLV) 
proves to be sensitive to adopted opacities and 3D radi¬ 
ation transfer capabilities of the spectral synthesis pro¬ 
cedure, and corresponding observations at inclinations 
from disc centre to close to the limb are very well fit by 
our 3D models (Secs. [3] and 0}. The good matching of 
absolute continuum intensity at disc centre (Sec. [3J con¬ 
firms that the effective temperature of the deep-lying 
continuum-forming layers is realistic, but does not pro¬ 
vide tight constraints for higher layers. The CLV of con¬ 
tinuum intensity is a better tool for that purpose, and 
our results provide a quite satisfactory match (Fig [3]) of 
solar limb darkening data, with the implication that the 
theoretical average temperature stratification should be 
relatively close to the correct one. The synthetic CLV 
values tend to be lower than the observations, in partic¬ 
ular at intermediate angles (/r z ~ 0.5), when considering 
intensities normalized to the values at disc centre, as is 
customary in CLV studies. So, this deviation adds up to 
the error in the absolute matching found in Sec. |3] (aver¬ 
age of 4%, see Fig. O. For an even more accurate fit of 
solar continuum data, lower opacities and/or increased 
radiative flux at ~ 1400 — 1800 nm (possibly due to us¬ 
ing a larger number of opacity bins in the simulations, 
or to atomic data improvements at those wavelengths), a 
mildly warmer temperature structure especially in inter¬ 
mediate to deep layers of the atmospheric model (thus, 
with a steeper temperature gradient at logTsoo nm — 1) 
and/or more magnetic flux, may be required. Note that 
precisely in the wavelength range in question the lin e-of- 
sight penetration in the Sun reaches its maximum ( Cray 
Il992l . chapter 9), which would indicate the need for bet¬ 
ter modelling in those deep, warm layers. Alternatively, 
the larger discrepancy present especially around 1500 nm 
(see Fig.[3l lower panel) may imply that improved obser¬ 
vations are required in that wavelength region. 

Regarding the [O 1 ] spectral lines targetted in this 
study, we obtain a good match of the synthetic data 
with the corresponding observations. In particular, when 
employing the 3D MHD “200 G” photospheric model, 
the match is achieved for logeMHD(C) 0 = 8.50 dex and 
log cmhd(O)© = 8.73 dex, while neither the HD nor 
the MHD models are able to fit observations if adopt¬ 
ing the lower solar abundances of lAsphmd et al.l ( 2009 1. 
loge(C) 0 = 8.43 dex and loge(O) 0 = 8.69 dex. To ap¬ 
preciate the effect of the magnetic fields with the largest 
possible clarity, in Sec. [5] we only showed results for the 
extreme cases of zero magnetic field and of the “200 G” 


MHD time-series, the latter being the case with the high¬ 
est vertical magnetic flux among our ensemble of simula¬ 
tions. Comparing these two cases, there is a clear effect 
of the modified average atmospheric temperature stratifi¬ 
cation on the line profile. However, the other MHD series 
(“50 G” and “100 G”) are more likely to be representa¬ 
tive of standard quiet Sun conditions. Additionally, for 
maximum accuracy when aiming at absolute solar chem¬ 
ical composition determinations, line broadening by the 
Zeeman effect should be included. Consequently, the fi¬ 
nal solar carbon and oxygen abundances that we could 
deduce from the present study may be somewhat differ¬ 
ent than obtained in Sec. [5] This needs to be checked 
in detail using a larger number of spectral lines - and, 
ideally, even better observational data - for the chemi¬ 
cal elements in question, in future studies. In any case, 
it is a significant result that, due to the sensitivity of 
the profiles to the temperature stratification, differen¬ 
tial effects are present for the two oxygen lines studied, 
and that with an appropriate choice of input chemical 
abundances, our MHD models can provide a good match 
of the corresponding observed spectra. Magnetic effects 
are therefore non-negligible, and given the observational 
evidence of the presence of magnetic fields all over the 
Sun, they should be taken into consideration. One may 
achieve a match of the HD model to the observed line 
spectra b y rais ing the abundances a little from the level of 
lAsplund et al.l (|2009tl . not so much as needed to achieve 
a good fit when using an MHD case. For each particular 
case (HD or MHD) one should test if one could find a 
comparable level of consistency in the set of derived C 
and O solar abundance values as we have found here for 
the “200 G” MHD case. 

The fair match of different observational quantities is 
an important result both in terms of confirming the qual¬ 
ity of the MHD models we employ, as well as in relation 
to the conclusion presented in Paper [Hi] where its authors 
argued that it is still premature to consider a revision of 
solar abundances based on 3D MHD modelling. Interest¬ 
ingly, although our 3D magnetoconvection calculations 
employ the same code as used by them, our resulting 
photospheric models have a different temperature strati¬ 
fication than theirs. Differences in the number of opacity 
bins used, equation-of-state treatment, spatial resolution 
and temporal variation of the snapshots may contribute 
to those differences. Our models predict warmer temper¬ 
atures in the middle photosphere, and as a consequence, 
a less steep temperature gradient in layers where many 
of the spectral lines that play a role in solar chemical 
abundance determination form. This can have a large 
effect on the emerging radiation spectrum, given that 
the excitation or ionisation processes involved in spec¬ 
tral line formation are predominantly dependent on that 
gradient. The results of this paper, namely the notice¬ 
able changes in the line spectra when going from HD to 
MHD and the excellent match that can be achieved in 
continuum and lines for the MHD models, together with 
further result s in the literature like th ose of Papers H and 
H3 or those of lBeck et al.l (1 20131 . 1 20151 ). which indicate an 
accurate average temperature stratification for our mod¬ 
els with ( B z } between ~ 50 G and ~ 100 G, confirm 
the need to consider magnetoconvection processes when 
solving problems that have a sensitive dependence on the 
shape of spectral features. 
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7. CONCLUSIONS 

In this work, we have studied the lit of various spectral 
features obtained through our 3D MHD numerical mod¬ 
els of convection to observational constraints, finding, in 
general a very good agreement. This is particularly the 
case for absolute continuum intensity matching as well 
as for centre-to-limb variation in the wavelength range 
600-2400 nm. The adopted opacities and the 3D a poste¬ 
riori radiation transfer processing prove to be important 
for a good match of observed solar continuum intensity 
from disc centre to limb. We have also shown that it is 
possible to achieve a good fit of the difficult [O i] 557 nm 
and 630 nm spectral features in spite of the blends and 
pecu liarit ies o f t hese lines. A convection model with the 
lAsplund et al.l (120091) abundances [loge(C) Q = 8.43 dex 
and loge(0)g = 8.69 dex] does not yield a satisfactory 
match, neither for the HD nor for any of the MHD cases. 
A better fit can be obtained by modifying the input abun¬ 
dances; for instance, a very good fit is obtained for our 
MHD “20 0 G” mo_del usi ng abu ndance values higher than 
those of lAsplund et al.l (|2009f ) by several centidex. Also 
importantly, the presence of magnetic fields in the con¬ 
vection models is seen to cause non-negligible changes in 
the line profiles (especially in the [O i] 630-nm case). This 
effect occurs indirectly via the modification in the aver¬ 
age temperature stratification of the atmosphere in the 
line-forming layers. The line asymmetries are also seen 
to be marginally different in the HD and MHD cases, so 
we cannot exclude the possibility that the differences in 
the average flow fields between the HD and MHD cal¬ 
culations play a role as well in modifying the line pro¬ 
file shapes. Adding to the results in this paper those 
from previous publications on abundances (Paper IH: HH). 
or those dealing with thermodynamic propert ies o r po¬ 
larization s ignals from our numerical models (jBeck et al.1 
I2013l.l20l5lf . we conclude that one should not neglect the 
magnetic field of the convection cells even in quiet Sun 
regions when aiming at accurate spectral modeling and 
chemical abundance determination. 


and of the first release (2014-02-05) of the “VALD3” 
data, as available from the web interface of the 
Vienna Atomic Line Database (VALD) project at 
http: / / vald. astro. uni vie. ac. at /~'vald3/php / vald. php 
Finally, we thank the anonymous referee for the helpful 
comments. 
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